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ABSTRACT 

We analyze the various interstellar components of the Hn region Sh2-132. The 
main stellar source is the double binary system that includes the Wolf-Rayet star 
WR153ab. We use radio continuum images at 408 and 1420 MHz, and Hi 21cm line 
data taken from the Canadian Galactic Plane Survey, molecular observations of the 
12 CO(1-0) line at 115 GHz from the Five College Radio Astronomy Observatory, and 
available mid and far IR observations obtained with the MSX and IRAS satellites, 
respectively. 

Sh2-132 is composed of two shells showing radio continuum counterparts at both 
frequencies. The emission is thermal in nature. The estimated rms electron density and 
ionized mass of the nebula are n e ~ 20 cm -3 and Mhii — 1500M Q . The distribution 
of the CO emission shows molecular gas bordering the ionized nebula and interacting 
with it. The velocities of the molecular gas is in the range -38 to -53 km s — 1 , similar 
to the velocity of the ionized gas. 

The emission at 8.3 /xm reveals a ring like feature of about 15' that encircles the 
bright optical regions. This emission is due to the PAHs and marks the location of 
photodissociation regions. 

The gas distribution in the environs of Sh2-132 can be explained in a scenario 
where the massive stars in the region photodissociated, ionized, and swept-up the 
dense molecular material from the parental cloud through their strong stellar winds 
and intense UV photon flux. 

Key words: ISM: bubbles - ISM: molecules - ISM: Hn regions - ISM: individual: 
Sh2-132 - stars: Wolf-Rayet - stars: individual: WR153ab 



1 INTRODUCTION 

The interaction between massive stars and the interstellar 
medium leads to the formation of Hn regions, interstellar 
bubbles, and new stars. More specifically, the interaction 
between the molecular medium and the UV stellar radia- 
tion is at the origin of photodissociation regions (PDRs). 
Molecules are dissociated by UV photons in the range 6 eV 
to 13.6 eV. The gas component ranges from ionized atoms 
at the PDR surface to molecules deeper into the molecular 
region as the FUV flux is absorbed. The FUV flux illuminat- 
ing these regions range between 1^ Go ^ 10 6 , where Go is 
the intensity of the radiation field in the solar neighborhood 
(Habing 1968), while the Hi density lies between 10 cm -3 
< n H < 10 5 cm" 3 (Tielens & Hollembach 1985, Sternberg & 
Dalgarno 1989, Sternberg & Dalgarno 1995, Bakes & Tielens 
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1998). PDRs are formed at the edge of Hn regions, plane- 
tary nebulae, reflexion nebulae, and stellar forming regions. 
Forbidden spectral lines in the far infrared as well as molec- 
ular lines are observed towards PDRs (e.g. Hogerheijde et 
al. 1995, Kaufman et al. 1999). 

An important component of PDRs is interstellar dust. 
Its presence plays a major role absorbing stellar radiation 
and reemitting it in the mid and far infrared. Polycyclic 
aromatic hydrocarbons (PAHs), which emit at 3.3, 6.2, 7.7, 
8.6, and 11.3 /jm, are tracers of PDRs (e.g. Watson et al. 
2008). 

In this paper, we study the interstellar environment 
of the Hn region Sh2-132 (Sharpless 1959), located in the 
Perseus spiral arm, paying particular attention to the pres- 
ence of PDRs at the interfase between the ionized and molec- 
ular gas. 

Sh2-132 was detected in Ha, [O III], and [S II] lines by 
Heckathorn et al. (1982). The nebula is centered approxi- 
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mately at (Z,6) ~ (102°50', -0°42') and presents two bright 
regions, which are surrounded by diffuse emission regions. 
They can be identified in Fig. 1, which displays the DSS-R 
image of the nebula. 

Churchwell & Walmsley (1973) performed radio con- 
tinuum observations at 2695 MHz with an angular resolu- 
tion of 18.'2, and, using measurements at several frequen- 
cies, derived a thermal spectral index. Harten et al. (1978) 
investigated the characteristics and morphology of the neb- 
ula based on radio continuum observations at 610 MHz us- 
ing the WSRT with a synthesized beam of 60". Most of 
the emission at this frequency originates from two optically 
bright regions named Shell A and Shell B, located at (I, b) ~ 
(102°57', -0°42') and (l,b) ~ (102° 47', -0°42'). Both shells 
coincide with the bright optical regions within the nebula 
and are indicated in Fig. 1. 

According to Harten et al. (1978), the optically bright 
section of Sh2-132 (at I > 102° 43') is an ordinary Hn region 
excited by the massive O-type star BD+55°2722 and the 
WR star WR153ab (= HD 211853), which are located close 
to the center of the nebula and linked to Shell B. Chu et al. 
(1983) reached to similar conclusions, based on the presence 
of the stars near the center of the nebula and on the absence 
of velocity changes in the nebular spectra. The main param- 
eters of these stars, as well as of other OB stars listed by 
Harten et al. (1978) that appear projected onto the Hn re- 
gion, are summarized in Table 1, which lists their galactic co- 
ordinates, spectral classification, visual magnitude V, colour 
indices (B — V) and (B — V)o, visual absorption Ay, absolute 
magnitude My, and derived spectrophotometric distance d. 
Uncertainties in the derived spectrophotometric distances 
of some of the stars were estimated using My values from 
Vacca et al. (1996) and Landolt-Bornstein (1982). The large 
distace derive for LS III+55°37 casts doubts on their relation 
to the nebula. Note that, based on photometric data, Panov 
& Seggewiss (1990) concluded that HD 211853 is a system 
with two pairs of stars, both hosting a WR component and 
an O-type star. 

Chu & Treffers (1981), based on Fabry Perot observa- 
tions of the Ha line, found that the ionized gas has LSR 
velocities in the range -53 to -45 km s~ , in agreement with 
previous results from Georgelin & Georgelin (1976, -50.9 km 
s" 1 ), Reynolds (1988, -48±1 km s" 1 ), and Fich et al. (1990, 
-47.1 km s _1 ). Quireza et al. (2006) detected ionized gas 
at -50.46 and -49.48 km s _1 from He and H radio recom- 
bination lines, respectively. Harten et al. (1978) estimated 
an emission measure EM = 3000 pc cm -6 for the entire 
nebula. 

New Ha and [O III] images taken by Miller & Chu 
(1993) revealed an arc of 4' in size to the south of the 
WR star (at [l,b] = [102°46', -0°43']). They concluded that 
the filament probably indicates the position of an ionization 
front and considered it to be a ring nebula related to the 
star. Based on optical data, Esteban & Rosado (1995) found 
velocities in the range -41 to -59 km s _1 for the arc and de- 
rived an electron density of 290±100 cm~ 3 from line ratios. 
The [OIII]/Ha ratio is consistent with values derived for ring 
nebulae, suggesting that WR 153ab is the main source of ion- 
ization (Esteban & Rosado 1995). These authors classify the 
arc as R s -type based on the morphology of this feature and 
on the absence of a clear evidence of expanding motions in 
the line profiles. This last point suggests that the massive 
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Figure 1. DSS-R image of the brightest section of Sh2-132. The 
grey scale is arbitrary. The large cross marks the position of the 
WR star, while smaller crosses indicate the location of the OB 
stars. 

progenitor of the WR star contributed in the shaping of the 
nebula (see Chu et al. 1981). 

Distance estimates for the WR star based on different 
methods range from 2.75 to 5.0 kpc. Massey (1981) de- 
rived a spectrophotometric distance of 5.0 kpc, while van 
der Hucht (2001) estimated 2.75 kpc based on its associa- 
tion to Cep OBI, with a distance uncertainty of 40%. Foster 
& Routledge (2003), using a method based on Hi column 
densities, estimated 3.2±0.5 kpc. A distance of 3.68 kpc was 
adopted by Churchwell & Wamsley (1973), based on the 
WN6+BOI spectral type of the WR star, known at that 
time. Georgelin & Georgelin (1976) derived a distance of 3.6 
kpc, while Fich, Blitz & Stark (1989) obtained 4.2±1.5 kpc. 

Adopting a mean LSR velocity of -48 km s _1 for the 
ionized gas, and taking into account the presence of non- 
circular motions in this section of the Galaxy (see Fig. 2a 
by Brand & Blitz 1993), a kinematical distance d k = 3.5±1.0 
kpc can be predicted, in close agreement with distances de- 
rived optically for the WR star. We adopt this distance in 
the present study. 

With the aim of investigating in some detail the nature 
of Sh2-132, we analyzed the distribution of the different gas 
components in its environs. To accomplish this goal we per- 
formed a multi-frequency study using radio continuum data 
at several frequencies, Hi 21cm line data, 12 CO(1-0) line 
data, and IR observations at different wavelengths. These 
data revealed that photodissociation regions are present at 
the interface between the ionized and molecular gas. 



2 DATABASE 

Neutral hydrogen 21cm line data from the Canadian Galac- 
tic Plane Survey (CGPS, Taylor et al. 2003) was used to 
analyze the Hi gas distribution, while 12 CO(1-0) data at 115 
GHz from the Five College Radio Astronomy Observatory, 
allowed us to study the molecular distribution. Intensities 
are shown as main beam brightness temperatures T m b- The 
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Table 1. O-B type stars towards the bright section of Sh2-132 



Star 




(1,6) 


ST 


V 


(B-V) 


(B - V) 




My 


d 












mag 


mag 


mag 


mag 


mag 


kpc 


HD 211853 


102 


3 46.'8 


-0°39' 


WN6/WC+O6IM 


9.08( 2 ' 


0.27< 2 ) 


(2) 


2.28( 2 ) 


-6.4( 2 ) 


2.75( 2 ) ,5.0W 


(WR 153ab) 






















LS III+55°37 


102° 


46f8,- 


-0°40;2 


07V< 3 ' 


11.01 


0.32(3) 


-0.29( 6 ' 


1.86 


-5.2( 5 ) 


7.5 


BD+55°2722 


102° 


48.'6,- 


-0°40;2 


08.5V< 3 ) 


9.91' 3 ) 


0.45(3) 


-0.28( 6 ) 


2.26 


-4.7(5) 


3.0 


LS III+55°40 


102° 


56f4,- 


-0°A2'.6 


09V( 3 ) 


12.08( 3 ) 


0.64(3) 


-0.28( 6 ) 


2.85 


-4.5(5) 


5.6 


LS III+55°42 


102° 


56f4,- 


-0°43;8 


B (7) 


12.38 


0.52 


-0.26( 6 ' 


2.42 


-2.45, -5.l(* 5 > 


5-0-10.3 


LS 111+55° 39 


102° 


57.'6,- 


-0°38f4 


B< 7 ) 


11.74 


0.45 


-0.26( 6 ) 


2.20 


-2.45, -5.l(* 5 ) 


2.5-8.5 


LS III+55°45 


102° 


59.'4,- 


-0°45?6 


08( 3 ) 


10.45(3) 


0.70(3) 


-0.27( 5 > 


3.01 


-4.9, -5.8(** 5 ) 


2.9-4.9 


LS III+55°49 


103° 


02.'4,- 


-0°48?6 


B< 7 ) 


11.50 


0.57 


-0.26( 6 ) 


2.5 7 


-2.45, -5.l(* 5 > 


1.9-6.4 



References: (1) Smith et al. (1998), (2) van der Hucht (2001), (3) GOS Catalog V2.2, Sota et al. (2008), (4) Massey (1981), (5) 
Landolt-Rormstein (1996), (6) Wegner 1994, (7) SIMBAD database, (*) Spectral classsification B0III-B2V was adopted, (**) spectral 

class V and III vas adopted. 



Table 2. Observational parameters of the Hi and 12 CO data 



00 00 FT 





Hi 


12 CO(1-0) 


Original FWHM 


I.'18x0.'98 


~1' 


Final FWHM 


l'Axl'A 




Vel. range (LSR) (km s" 1 ) 


-165.+57 


-165, +57 


Vel. resolution (km s _1 ) 


1.32 


0.824 


RMS noise T b (K) 


1.7 


0.15 



Hi data were smoothed to l'A x l'A to facilitate the identifi- 
cation of the Hi features. The main observational parameters 
of these data bases are detailed in Table [2] 

Radio continuum observations at 408 and 1420 MHz 
were also extracted from the CGPS. The image at 408 MHz 
was obtained with a synthesized beam and an rms noise of 
3l4x2.'8 and ~1.1 K, respectively; the corresponding values 
for 1420 MHz are l.'0x0!81, and 0.063 K. 

Mid-infrared (MIR) images were extracted from the 
MSX Galactic Plane Survey (Price et al. 2001) at 8.3/xm 
(A-band), 12.1/xm (C-band), 14.7rtm (D-band), and 21.3/xm 
(E-band), with an angular resolution of 18'.' 4. MIR and far- 
infrared (FIR) images at 12, 25, 60, and 100 /im from the 
IRAS satelite (HIRES data) are also included. These images 
have angular resolutions in the range 0.'5 to 2!0. 



3 RESULTS 

3.1 Radio continuum emission 

Figures [5] and show the CGPS radio continuum emission 
at 1420 MHz. The grayscale in Fig. 2 has been chosen to em- 
phasize the large-scale region of faint emissivity to the north 
of WR 152 (Cappa et al. 2010, hence the "burnt-out" ap- 
pearance of the region around WR 153ab). The emission dis- 
tribution at 1420 MHz resembles an upside-down "U". This 
semicircular structure consists of two areas of strong emis- 
sion coincident with Shells A and B detected by Harten et al. 
(1978), centered at (I, b) ~ (102°34'42'.'4, -00°46'9'.'8) and at 
(l,b) ~ (102°26'35'.'8,-00°43'10'.'2), respectively. WR153ab, 
indicated by the larger cross in Fig. 1, appears projected 
onto the central region of Shell B. Thus, the bright radio 
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Figure 2. Full resolution CGPS image at 1420 MHz emphasizing 
faint emission to the north of WR 152 (indicated by a cross 
symbol in the lower right part of the figure). The arrow points 
to a continuum feature corresponding to the large HI structure. 



emission at I ^ 102? 5 is in excellent agreement with the op- 
tical emission (see Fig. [T]) and with the image at 610 MHz by 
Harten et al. (1978). Both Shells A and B can be identified 
at 1420 MHz. Note in particular the small arc of emission at 
(l,b) = (102?77, -0?73) (best identified in Fig. 3), which is 
the radio counterpart of the 4'-arc detected at optical wave- 
lengths by Miller & Chu (1993). 

The faint structure centered at (I, b) = (102?2, -0?83) is 
related to WR 152 and has been analyzed in a previous pa- 
per (Cappa et al. 2010). Of particular interest is the curved 
filament indicated by an arrow in Fig. 2. This filament can 
be followed from b » -0?75 to nearly b w -0?05. The fig- 
ure also shows a faint barely visible northeastern section. 
This structure is likely the radio continuum counterpart of 
the ~60' diameter Hi structure discussed in Cappa et al. 
(2010). 

Three bright point sources, best seen on Figure [3] are 
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Figure 3. Full resolution 1420 MHz CGPS image showing the 
eastern and western regions of Sh2-132 used for TT-plot analyses. 
The radio sources in direction to Sh2-132 are indicated by 10, 13 
and 14. The white contour line marks the 9 K contour of the 
convolved 1420 MHz image used for the analysis. See text for 
details. 



projected onto the nebula near Shell A. Their coordinates, 
flux densities and identifications are included in Table 3. 
The spectral index estimates indicate that sources 10 and 
14 are non-thermal in nature and probably extragalactic, 
while source 13 is thermal. The value of the spectral index 
indicates that the radio source is optically thick at least at 
408 MHz. 

The CGPS image at 408 MHz does not show as many 
details as the one at 1420 MHz, as it is of much lower res- 
olution. However it can be used to study the spectral in- 
dex distribution of the radio continuum emission around 
WR 153ab. As a first step, we removed the point sources 
which are present over the region of interest and convolved 
the 1420 and 408 MHz images to the same spatial resolution 
of 3'A. We then constructed a TT-plot, in which the bright- 
ness temperature T), at one frequency is plotted point-by- 
point against the brightness temperature at the other fre- 
quency. The brightness temperature spectral index /3, where 
Tf, oc is 13 , is directly related to the slope of the straight line 
fitted by regression, /3 = - log(slope)/ log(1420/408). The 
usual flux density spectral index a (S v oc u a ) is simply 
a = P + 2. The TT-plot method is immune from errors 
in the zero-level of the individual images (which the sim- 
pler technique of dividing one map by the other, even after 
background subtraction, is not). 

Three different TT-plots were constructed, one for the 
entire Sh2-132 region and one for the eastern and west- 
ern halves, which correspond to the two shells identified by 
Harten et al. (1978). Figure[3] which is the full resolution im- 
age at 1420 MHz, shows these different regions, the dashed 
lines delimiting the eastern and western parts. The three 
bright small diameter sources coincide with sources 10, 13 
and 14 in Table 2 of Harten et al. (1978) and had been 
subtracted from the original image before convolution to a 
circular beam of 3'A. 

The results of the TT-plot analysis are presented in 



a = - 0.03 ± 0.03 




Figure 4. TT-plots corresponding to (a) the whole structure, (b) 
the eastern half and (c) the western half. For all plots, a back- 
ground of 9 K and 80 K has been subtracted from the convolved 
images at 1420 and 408 MHz, respectively. 
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Figurc[4] In each case, the analysis was limited to regions 
brighter than 9 K at 1420 MHz and 80 K at 408 MHz. Al- 
though a significant amount of scatter is present on the dia- 
gram, especially for the eastern sector, there is no evidence 
for the presence of significant variations of the spectral index 
and a further study of spatial variations is not warranted. 
The best-fit spectral index is a = -0.03±0.03, corresponding 
to thermal emission. 

In order to obtain some radio continuum parameters, we 
first derive the total number N u of UV ionizing photons per 
second required to produce the observed free-free emission. 
This can be written as (e.g. Chaisson 1976): 



(2 



N u = 0.76 x 10 47 T 4 -°' 45 ^hz S v dl pc 



(1) 



where S u is the flux density in Jy, T4 is the electron temper- 
ature in units of 10 4 K and i/ghz the frequency in GHz. We 
use the CGPS image at 1420 MHz. The derived flux den- 
sity is significantly dependent on how the source dimension 
is defined. The surrounding background is somewhat vari- 
able, nevertheless it is reasonable to assume that the source 
emission corresponds to brightness temperatures above 6.5 
to 7 K. Using an average background of 7 K, we obtain 
a flux density at 1420 MHz of 23 ± 5 Jy, the error being 
determined from a second estimate using a background of 
6.5 K. From the above equation for N u and taking T4 = 0.8 
(Quireza et al. 2006), we obtain N u = 2.2 x 10 49 s" 1 . Consid- 
ering two 061-type stars and two WN6 as components of the 
multiple system, and the spectral type from BD+55°2722, 
LSIII+55°45, and LS+55°40 (see Table 1, related to the Hll 
region, we obtain the UV ionizing photons per second emited 
by the stars, Qo=(16.6±3.0) x 10 49 s _1 , for solar metalicity 
from Smith et al. (2002). 

The comparison between the number of UV photons 
emitted by the stars Qo and the number of UV photons 
used to ionize the gas N u indicates that the massive stars in 
the nebula are responsible for the ionization. A large frac- 
tion of the UV photons are also used to heat the associated 
interstellar dust. This result differs from the one by Harten 
et al. (1978) since their estimate of Qo was too low. 

We can now use the models of Mezger & Henderson 
(1967) to infer the properties of the extended ionized gas, 
namely its mass Mhii and the rms electron number density 
n e . Letting 0Ga be the observed half-power width of the 
extended emission (in minutes of arc, after correction for 
the beam), we obtain: 

n e /cm- 3 = 6.35X 10 2 m T 4 175 v° G £ S° v 5 9^ d^{2) 

M HII /M = 0.386 u 2 a~ 1/2 T 9 ' 175 u°o°l S 9 ' 5 6^1 d 2 kp \ (3) 

where a(u,T e ) is the Gaunt factor (about unity for our pur- 
poses), dkpc is the distance in kpc and the parameters ui and 
U2, both of order unity, depend on the assumed model for the 
source. These parameters are given by Mezger & Henderson 
(1967) for three different idealized models. From Figur^U it 
is apparent that the nebulosity can be roughly broken up 
into two different regions separated by a nearly vertical line 
running at a longitude of ~102?87. These two regions cor- 
respond to shells A and B from Harten et al. (1978). The 
flux densities of the two sub-regions are 11 ±2 Jy and 12 ±3 
Jy, respectively. For simplicity, we model each half-region as 
a sphere of constant density and of angular size 12'. From 
Mezger & Henderson (1967), we find ui = 0.78, u 2 = 1.29. 
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Figure 5. Upper pannel Sh2-132 at 8.3 /im (MSX band A). The 
grayscale goes from 8.0xl0 -7 to 3.5xl0~ 6 W m~ 2 ster -1 . Ex- 
tended filaments and point-like sources are indicated (see text). 
Bottom pannel: Overlay of the DSS-R image and the image at 
60/im. IR contours are 50 to 190 MJy ster -1 in steps of 20 MJy 
ster -1 , and from 300 to 900 MJy ster -1 in steps of 200 MJy 



Assuming T4 = 0.8, from which we determine a -1 / 2 = 1.04, 
we finally have an rms n e = 20 cm -3 and Mhii = 15OOM0 
for each region, where, given the error in estimating fluxes, 
we have simply taken S» = 10 Jy for each half. The derived 
rms electron density is compatible with estimates by Harten 
et al. (1978) for Shells A and B. 

3.2 IRAS and MSX emission 

The IRAS (60/xm) and MSX band A (8.3/xm) emissions are 
shown in Fig[5] Although we analyzed the four IRAS bands, 
we include the 60 pm-image only since the extended emission 
is clearly seen at this wavelength. Emission at 8.3/im (up- 
per panel) in diffuse nebulae has a contribution from both 
PAH and dust continuum. As pointed out by Povich et al. 
(2007), in very energetic regions like M 17, dust continuum 
can be very strong. However, for the case of Sh2-132, the as- 
sumption that MSX band A is dominated by PAH emission 
seems to be reasonable. Consequently, emisison in this band 
indicate the presence of photodissociation regions, while the 
emission in the far IR originates in large dust grains. 
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Table 3. Small diameter radio sources towards Sli2-132. 



# 


I 


Sl420 


S408 


a 


Sl78 


S365 


S610 


S1420 


S4850 


Identification 




b 


mJy 


mjy 




mJy 


mJy 


mJy 


mJy 






(1) 




(2) 


(2) 


(2) 


(3) 


(4) 


(5) 


(6) 


(5) 




10 


102°59'50'.'90 


195±10 


550 


-0.8 




819± 87 


330 




<40 


7C 2218+558( 7 > 




-0° 40' 23'.' 35 




















13 


103°3'30'.'24 


365+15 


280 


+0.2 






160 


222.6±7.5 


<40 


NVSS J222034+561438 




-0°41'28'/14 




















11 


103°4'40'.'04 


450±20 


1200 


-0.8 


3100±500 


1738±62 


810 




60 


TXS 2219+559 




-0° 47' 27? 04 





















References: (1) Number from the list of Harten et al. (1978); (2) this paper, derived using CGPS data; (3) Gower et al. (1967); (4) 
Douglas et al. (1996); (5) Harten et al. (1978); (6) Condon et al. (1998); (7) Waldram et al. (1996). 



The image at 8.3 /im (upper panel) reveals the existence 
of a ring-like feature centered at (l,b) = (102°53', -00°40') 
of about ~15' in size. This feature can be also identified at 
60 /im (lower panel). The comparison of the optical emission 
with the emission at 8.3 /mi shows that the mid IR emission 
encircles most of the nebula, following the outer rim of the 
ionized gas. This fact is more noticeable for Shell B and is 
particularly striking near (l,b) = (102° 47' -00° 43'), where 
the 4' ring nebula is present. The bottom panel shows that 
the emission distribution at 60 /mi resembles that in the 
optical. 

In order to estimate some dust parameters, we have di- 
vided this ring-like structure in three main filaments, named 
fl, f2, and f3, which are indicated in the upper panel of Fig. 
5. Note that f3 borders the region with high extinction seen 
in the optical image near (l,b) = ( 103° 0' -00° 47'). 

A number of interesting point-like sources located onto 
the filaments can also be identified in the MSX and IRAS 
images. They are indicated in the figure as psl, ps2, and 
ps3. Infrared source psl at (l,b) = (102°39'45",-00°36'24") 
is the MSX source G102. 6627-00. 6067, coincident with 
IRAS 22160+5551; ps2 at (l,b) = (102°48'18",-00°43'6") 
corresponds to G102. 8051-00. 7184 and IRAS 22172+5549; 
while ps3 at (l,b) = (103°03'44",-00°41'27") is G103.0624- 
00.6911 and IRAS 22187+5559. 

After background substraction, the IR flux densities 
of the filaments in the four IRAS bands and in the MSX 
band A were calculated. These values are listed in Table 2] 
The quoted errors reflect the uncertainties in the determi- 
nation of the background. Emission in bands C, D, and E is 
very weak in comparison with the strong background emis- 
sion leading to large uncertainties in flux density estimates. 
Values for these bands are not included. The three point 
sources, whose flux densities in the MSX and IRAS bands 
are included in the table, can be classified as YSO candidates 
following criteria by Junkes et al. (1992) and Lumsden et al. 
(2002) for IRAS and MSX point sources, respectively. In the 
case of the IRAS criteria, Junkes et al. (1992) consider (i) 
S100 >20 Jy, (ii) 1.2< Swo/S 6 o <6, (iii) S 60 /S 25 >1, and 
Q6O+Q100 Js4, where Si is the flux density in the bands cen- 
tered at 25, 60 and 100 /im, and Qi is the quality of the 
fluxes at these bands. Lumsden et al.'s (2002) criteria are 
(i) S21/SS > 2 and (ii) S14/S12 > 1, where Si is the flux 
density centered at 8.3, 12.1, 14.7, and 21.3 /im. 

It has been demonstrated (Scoville & Kwan 1976) that 
the infrared spectrum of radiatively excited dust grains can 
be fitted by a modified Planck function B(Td)iy m , were v m is 



the grain emissivity and m increases from 1 to 2 depending 
on the wavelength and composition of dust grains (Schwartz 
1982). Adopting m = 1.5, according with the molecular envi- 
ronment (Whittet, 1992), and using the measured flux den- 
sities at 60 and 100 /im, we derived the dust color tempera- 
ture Td summarized in Table(4] IR luminosities for the point 
sources are also included. Infrared filaments are excited by 
the intense stellar radiation field of WR 153ab and the other 
massive stars in the region. The slightly higher dust color 
temperature derived for fl can be explained bearing in mind 
that this filament is closer to the exciting stars than f2 and 
f3. In Fig. [6l the MIR and FIR spectral energy distribu- 
tion corresponding to the filaments are shown. Note that the 
measured fluxes at shorter wavelengths exceed the values of 
our model, possibly as the result of the emission of nebular 
lines which may contribute to the infrared flux densities. 



3.3 Molecular gas distribution 

The analysis of the 12 CO(1-0) emission distribution leads to 
the detection of two structures with velocities in the ranges 
from -53.4 to 46.0 km s"\ and from -46.0 to -38.6 km s" 1 . 
Fig. shows the averaged CO emission (contours) for the 
two molecular structures superposed onto the emission at 
8.3 /im. 

The CO feature present between -46.0 and -38.6 km 
s _1 (upper panel) exhibits a ring-like appearance between 
I ~ 102° 35' and 102° 50', while for I > 102° 45', the emission 
is concentrated along b ~ -0°35'. The whole structure will 
be referred as cloud A. Two smaller molecular clouds are 
detected at (l,b) ~ ( 103° 2' -0° 30') (indicated in the figure 
as cloud C) and at (l,b) ~ ( 103° 4', -0° 42') (cloud B). The 
emission at 8.3 /im correlates very well with the molecular 
emission. In particular, the brighter emission regions in the 
mid infrared coincide with clumps in cloud A [for example at 
(l,b) ~ ( 102° 43' -0° 35'), ( 102° 44' -0° 43'), ( 102° 47' -0° 44'), 
and ( 103° 5' -0° 42')]. Cloud B coincides with the strong ra- 
dio source NVSS J222034+561438, which is an YSO can- 
didate. The comparison of the ionized and molecular gas 
distributions shows that the molecular emission encircles 
the brightest optical emission regions analyzed by Sanchez- 
Monge et al. (2008). 

The lower panel displays the molecular emission be- 
tween -53.4 to -46.0 km s" 1 . The CO emission is clumpy, 
being concentrated in clouds D, E and F. Cloud C is also 
detected in this velocity range. Cloud E overlaps a section of 
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Table 4. IR filaments and point sources. 



Source 








Fluxes [Jy] 








T d (K) 


Liras(Lq) 




8.3 /an 


12.1 fim 


14.7 (im 


21.3 fim 


12 (im 


25 fim 


60 fim 


100 fim 






fl 


46±1 








72±4 


192±24 


1370±130 


2765±145 


32±7 


43800 


f2 


16.0±3.5 








35±6 


77±11 


484±51 


1201±160 


29±6 


16900 


f3 


1.8±0.3 








4.9±0.5 


32±11 


133±40 


343±37 


32±6 


5000 


psl 


1.45 


1.85 


0.97 


2.14 


3.03 


4.28 


74.1 


274 


28±5 


2700 


ps2 


2.22 


3.93 


3.78 


7.02 


7.8 


16.1 


106 


318 


30±5 


3900 


ps3 


1.0 


1.6 


1.40 


4.26 


7.02 


17.2 
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36±8 


6800 
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Figure 6. IR spectral energy distribution for the IR filaments. 
The lines in the diagram corresponding to fl, f2, and /3 show the 
emission of a blackbody with a T = 32, 29 and 32 K, respec- 
tively, and a y 15 -emissivity law. The boxes denote flux densities 
in the four IRAS bands, and the circle, that corresponding to the 
MSX band A. 



cloud A and follows the ring-like feature identified at 8.3 fjxci. 
Cloud F coincides with part of the IR filament at (l,b) ~ 
(102°44' -0°42'). The section of Cloud D running from (/, b) 
= (102°56',-00°51') to (l,b) = (103°02' ,-00°44') coincides 
with f3. Cloud D also displays a clump at (I, b) ~ (103°05',- 
0°44'), which coincides with a region of high optical extin- 
cion and is probably connected with cloud B. 

Table [5] gives the physical parameters of each molecu- 
lar cloud. The masses of the different clouds were estimated 
by integrating the CO line intensity as W(CO) = J T mb dv 
= T m t,Av, where T m i is the average main beam brightness 
temperature of the molecular cloud over the velocity in- 
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Figure 7. Overlay of the averaged 12 CO emission distribution 
for the molecular clouds between -46.0 to -38.6 km s — 1 (upper 
panel) and -53.4 to -46.0 km s — 1 (lower panel) and the MSX 
band A emission (grayscale). Contours correspond to 3, 7, and 
15cr. The grayscale goes from 0.1 xlO -7 to lxlO -7 Watts m -2 
ster . The crosses indicate the positions of the massive stars. 



terval in which the cloud is observed. To estimate the H2 
column density N(H 2 ), the relation X = N(H 2 )/W(CO) 
= (2.3±1.2)xl0 20 mol cm" 2 KT 1 km" 1 s (Grenier & Le- 
brun 1990) was used. The molecular mass was derived from 
M(M ) = 4.2xl0 20 N(H 2 )d 2 Q, where d is the distance in 
pc and Q is the solid angle in steradians. Errors in T mb and 
in N(H 2 ) arise in background uncertainties, while errors in 
masses come from background and distance uncertainties. 
The adopted distance is 3.5±1.0 kpc. 

To determine the volume density, we assumed a spheri- 
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Table 5. Physical parameters of the molecular clouds. 









Clouds 








A 


B 


C 


D 


E 


F 


Average T mb (K) 


1.30±0.30 


0.90±0.15 


0.60±0.15 


1.20±0.30 


1.60±0.15 


0.70±0.15 


Velocity range (km s _1 ) 


-38.6,-46.0 


-38.6,-46.0 


-38.6,-53.4 


-46.0,-53.4 


-46.0,-53.4 


-46.0,-53.4 


N(H 2 ) column density (10 20 cm~ 2 ) 


6.50±1.60 


3.30±0.60 


4.20±1.05 


5.90±2.40 


9.0±0.90 


1.95±0.20 


H 2 mass (10 3 M Q ) 


4.25±1.90 


0.80±0.35 


0.04±0.02 


1.30±0.60 


1.55±0.70 


0.95±0.45 


Volume density (cm" 3 ) 


37±22 


2150±1290 


110±65 


95±55 


115±70 


320±190 


N(HI) column density (10 21 cm~ 2 ) 


3.7±0.9 


3.5±0.5 


3.7±0.9 


3.3±0.8 


3.4±0.3 


3.3±0.6 


Visual absorption A v (mag) 


2.5±0.8 


2.2±0.5 


2.8±0.7 


2.4±0.7 


2.8±0.6 


2.7±0.5 



cal geometry for all the clouds, except for Cloud A, for which 
an elongated geometry was adopted. 

Figure [8] shows two velocity-position maps of the CO 
emission. We have plotted Galactic longitude in the y-axis 
since the molecular morphology is mostly elongated in this 
direction. The map in the upper panel is the result of 
integrating within A6 = (-0°55'30" ,-0°39'54"), while the 
map in the bottom panel was integrated within A6 = (- 
0°39'54" ,-0°29'06"). The upper panel shows two clouds at 
I ~ 102°45' and I ~ 103°00', having velocities in the range 
Av ~ (-40,-45) km s _1 and (-45,-50) km s _1 , respectively 
(two small cloudlets are also detected at v ~ -30 km s _1 at 
I ~102 °50' and I ~103 °03'). Clouds F and D are clearly 
identified. The bottom panel shows cloud A extending from 
-102° 35' to 103° 13'. The velocity range goes from -40 km 
s _1 at I ~ 103° 10' to -50 km s _1 at I ~ 102°35'. The image 
shows some small cloudlets at I ~ 103° 23' at v ~ -38 km 
s _1 , I < 102°40' with velocities in the range ~ -8 to -5 km 
S _ , and at I > 103°25' at v ~ km s" 1 . Molecular gas 
with velocities higher than -10 km s _1 is unconnected to 
Sh2-132. 

The agreement in velocity between the ionized and 
molecular gas and the fact that the molecular emission bor- 
ders the brightest ionized regions indicate that molecular 
gas with velocities between -38 and -54 km s _1 is associ- 
ated with Sh2-132 and the massive stars in the region. 

We adopt -45 km s~ as the systemic velocity of the 
molecular gas. Taking into account the non-circular motions 
in this section of the Galaxy shown by Brand & Blitz (1993), 
a kinematical distance dk = 3.0±1.0 kpc can be estimated, 
compatible with the distance adopted in SectfJ] 

3.4 The distribution of the Hi emission 

The upper panel of Fig[9] displays the Hi emission distribu- 
tion in the velocity range where the CO emission is present, 
i.e. between -31.1 and -58.4 km s _1 . The white circle has 
the dimensions of the Hn region. Neither a region lacking Hi 
gas coincident in position with Sh2-132 as a whole nor an 
Hi shell surrounding it are evident in the image, although 
the small cavity near (l,b) ~ ( 102° 58' ,-0° 40') is probably 
linked to shell A. 

However, a gradient in the Hi emission distribution is 
evident for b < -0°10'. The lower panel displays the emission 
corresponding to I = 102°50' from b = 0°00' to -1°30' from 
the image in the upper panel. Tj, and N(HI) are shown on 
left and right y-axes, respectively. The image shows that 
the Hi brightness temperature is Tf, ~ 85 K at b ~ -0°15', 



while 65% this value is detected at 6 ~ -1°30'. The fitted 
curve, showed as a dashed line, confirms the presence of a 
IVgradient. 

The fitting curve was used to estimate the contribu- 
tion of the Hi gas to the visual absorption A v . We calcu- 
lated the visual absorption A v due to the molecular and 
neutral hydrogen related to the Hn region by using the rela- 
tion (N(HI) + 2N(H 2 ))/E(B-V) = 5.8 xlO 21 atoms cm~ 2 
mag -1 (Bohlin et al. 1978), where N(HI) and N(H 2 ) are the 
Hi and H2 column densities, respectively, while E(B — V) is 
the color excess. We adopted E(B — V) = A v /3.05 and read 
off N(HI)(b) from the diagram. The estimated Hi column 
densities, as well as the derived A v -values for each molecular 
cloud are included in Table [S] 



4 PHOTODISSOCIATED REGIONS 

The presence of PAH emission at the interfase between the 
ionized and molecular materia, e.i. at the border of many 
Hn regions, suggests that molecular gas is being photodis- 
sociated by the UV photons emitted by the massive stars. In 
the following subsections we analyze the presence of PDRs 
in two selected areas within the Hn region, one around 
WR 153ab and BD+55°2722, linked to shell B, and the other 
around LS+55°39, linked to shell A. 



4.1 Region near WR 153ab and BD+55°2722 

The upper panel of Figure 10 displays an overlay of the 
DSS-R image and the CO emission distribution for a small 
region of Sh2-132 near WR153ab and BD+55°2722. Thin 
contours correspond to Cloud A with velocities in the range 
-38.6 to -46.0 km s" 1 , while thick ones delineate Clouds E 
and F having velocities between -46.0 and -53.4 km s _1 . 
The crosses mark the position of the stars. 

The morphological agreement between the optical emis- 
sion and the inner borders of the molecular clouds is ex- 
cellent, particularly for the 4' optical filament at (I, b) ~ 
(102°46'30", -0°43'), where the contour corresponding to 
2.25 K closely follows the optical emission. Note also that 
difusse emission is also encircled by molecular emission. The 
existence of emission at 8.3 |im bordering the ionized region 
(see Fig. 7) strongly indicates the presence of a PDR. 

To confirm the existence of a PDR in this region we 
compared the distribution of the ionized and molecular gas, 
and that of PAHs. After convolving the DSSR and MSX 



The ISM linked to Sh2-132 9 



102 40 



50 



1 

u 103 00 
o 
3 
3 

w 



I 



I 



20 



30 



Ab=[-55'30",-39'54"] 



102 40 



50 



1 

u 103 00 
o 
3 
3 

10 




00 10230 
GALACTIC LONG. 



T b (K) 




N HI (cm- 2 ) 



b (arcmin) 

Figure 9. Upper panel: Average Hi emission distribution between 
—31.1 and —58.4 km s —1 . The grayscale goes from 50 to 93 K, while 
contours correspond to 60, 70, 80, and 90 K. The white circle has 
the dimensions of the Hn region Sh2-132. Bottom panel: Diagram 
showing Tj, vs. b from b= —90' to 0' for I = 102° 50' from the image 
above. The dashed line shows the best fit to the data. 



20 



30 



Ab=[-39'54",-29'06"]_ 



_L 



-10 -20 -30 40 
Kilo VEL0-LSR 



■50 -60 -70 



Figure 8. Upper panel: Averaged 12 CO velocity distribution for 
Ab = (-55'30",-39'54"). The grayscale goes from 0.15 to 3 K. 
Contours arc 0.45, 0.75 and 1.05 K. Bottom panel: Averaged 12 CO 
velocity distribution between Ab = (-39' 54", -29' 06"). Contours 
and grayscale are the same as in the upper panel. 



band A images to the CO beam (46"), we averaged the opti- 
cal, molecular and band A emissions in concentric rings from 
the position of the WR star within an angle of 300° from 
(/,&) ~ (102°52',-0°38') clockwise up to (l,b) ~ (102° 49',- 
0°33'). The result of this average is displayed in the bottom 
panel of the figure, where intensities in the different bands 
arc normalized to the maximum value of each band. The 
values in the abscisae indicates distance from the WR star. 



The diagram reveals a clear stratification, with the optical 
emission peaking near the stellar position, and the emission 
at 8.3 fjxa peaking 3' far from the star. The molecular emis- 
sion increases up to 6' from the star, where the maximum 
is present. The observed distribution is typical of molecu- 
lar clouds exposed to UV radiation, which photodissociates 
them (van der Werf et al. 1996). 

4.2 Region near LS+55°39 

The region around LS+55°39 is particularly interesting. The 
upper panel of Fig. 11 shows that the star is close to a bright 
optical filament of ~2' in size. The fact that the optical 
filament partially follows the border of cloud A indicates 
that the UV photons of the star are photodissociating and 
ionizing the dense gas. The high extinction region at (I, b) ~ 
(102°57',-0°40') coincides with a section of cloud D. 

With the aim of investigating the presence of a PDR 
near the optical filament, we plotted the normalized in- 
tensity of optical, MSX band A, and CO emissions as a 
function of the distance to LS+55°39. The plotted emission 
corresponds to the region indicated by the black arrow in 
the upper panel. As for the region around WR 153ab and 
BD+55°2722, the optical emission peaks closer to the ex- 
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Figure 10. Upper panel: Overlay of the DSS-R image and the CO 
emission in the region near WR153 and LS+55°2722, indicated 
as crosses. Thin lines delineate Cloud A and thick lines, clouds E 
and F. Bottom panel: Normalized intensities of the optical, MSX 
band A, and CO emissions in concentric rings as a function of the 
distance to the WR star. 
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Figure 11. Upper panel: Overlay of the DSS-R image and the 
CO emission in the region near LS+55°39, indicated by a cross. 
Thin lines delineate Clouds A and D. Bottom panel: Normalized 
intensities of the optical, MSX band A, and CO emissions in con- 
centric rings as a function of the distance to the B star. 



citation source than the emission in MSX band A and the 
molecular emission, revealing the presence of a PDR. 

Taking into account that LS+55°39 is a B0V star with 
logQ(-ff°) = 46.23 (Vacca et al. 1996), and assuming that 
it is immersed in an ionized interstellar medium with a 
mean electron density n e = 20 cm" 3 , we can estimate the 
Stromgren radius r s of the Hn region created by the star 
from — ^4 = r s (Mezger & Henderson 1967), where olb 
— 2xl0~ 13 cm^V 1 . The result is r s = 5.3 pc. Bearing in 
mind that the angular distance from the star to the optical 
filament is 2.2', or 2.25T0.7 pc at 3.5T1.0 kpc, we can con- 
clude that the B-type star is responsible for ionizing the gas 
in the region and contributing to the photodissociation the 
the molecular gas. 



5 SCENARIO AND CONCLUSIONS 

In this paper we have analyzed the distribution of the ion- 
ized, neutral atomic and molecular gas and that of the in- 
terstellar dust in Sh2-132. 

This study revealed the presence of a number of molec- 
ular clouds encircling the ionized nebula, with a mass of 
9000 Mq. The velocity of the molecular gas, which is in the 
range from -38.6 to -53.4 km s _1 coincides with the veloc- 
ities of the Ha line, and He and H radio recombination 
lines, strongly indicating that the ionized gas is associated 
with the nebula and interacting with the ionized gas. 

The ring-like appearance of the gas and dust distri- 
bution in the environs of the WR star suggests the ac- 
tion of stellar winds, that sweep up and compress the gas. 
The mechanical energy E w released by the multiple system 
WR 153ab into the ISM can be roughly estimated taking 
into account the two WR components of the binary system 
during the duration of the WR phase of a massive star, the 
previous main sequence phase of the current WR stars, and 
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the two O-star components. We adopt a mass loss rate log 
(M) = -6.64 and a terminal velocity = 1950 km s _1 , 
corresponding to a 08.5V star (Smith et al. 2002) as pa- 
rameters of the stellar wind of the O-type components. For 
the WR components we adopt a conservative mass loss rate 
M =5x 10~ 6 M yr -1 and a terminal velocity Vao = 2000 km 
s _1 . Adopting lifetimes of 3xl0 6 yr and 0.5xl0 6 yr for the 
O and WR phases of the stars, respectively (Conti & Vacca 
1990, Meynet & Maeder 2005), E v results to be ~ 7.2xl0 50 
erg. 

Considering the molecular clouds E, F and a section 
of cloud A, and the dust and ionized components, with an 
expansion velocity of 7.4 km s _1 , the kinetic energy of shell 
B is ~2.80xl0 48 erg. This result indicates that the multiple 
system can shape shell B. 

The emission at 8.3 fim appears concentrated in a ring- 
like structure bordering the brightest sections of the nebula, 
indicating that molecular material is being photodissociated 
at the interface between the ionized and molecular gas. This 
configuration resembles the IR bubbles analyzed by Wat- 
son et al. (2008). To confirm the presence of photodissoci- 
ation regions we have analyzed in some detail two regions 
within the complex. One in Shell B around WR 153ab and 
BD+55°2722, and the other in Shell A around LS+55°39. 
The relative positions of the ionized gas, PAHs, and molec- 
ular gas in these regions show the stratification typical in 
PDRs. Thus, the distribution of H2 plays a main role in 
shaping the Hn region. 

The emission distribution at 60 /xm shows the same 
ring-like structure, along with emission inside the ring. The 
presence of emission in the far IR inside the Hn region, which 
originates from large dust grains, is compatible with the ex- 
istence of PDRs since PAHs are destroyed inside Hn regions, 
while large grains are not. 

The number of UV photons emitted by the massive stars 
in the region is seven times higher than the UV photons 
used to ionize the gas, indicating that the massive stars are 
responsible for the ionization of the region. Clearly, a large 
fraction of the available UV photons are used in the heating 
of the interstellar gas and in the photodissociation of the 
dense material. 

Different studies of Hn regions such as those towards 
RCW79 (Zavagno et al. 2006), RCW82 (Pomares et al. 
2009), and Sh2-212 (Deharveng et al 2008) show how the 
interaction between the massive stars with the surrounding 
molecular material induces star forming regions through the 
"collect and collapse" process (Elmegreen and Lada 1977) 
in the borders of these Hn regions. For example, the paper 
on RCW79 (Zavagno et al. 2006) reveals the formation of 
YSOs inside pre-existing molecular cores located at the bor- 
der of the ionized region. By using the analytical model by 
Whitworth et al. (1994) to determine the time at which the 
fragmentation occurs and the size of the structure at that 
time, they conclude that the collect and collapse process 
may not be the only one at work. Sh2-212 (Deharveng et al. 
2008) and Sh2-104 (Deharveng et al. 2003) represent good 
examples of stellar formation in which the molecular clouds 
surround the Hn regions and the PAH emission is located at 
the interfase between the ionized and the molecular gases. In 
particular, Sh2-104 is encircled by a molecular ring showing 
four dense molecular clouds containing many dense cores. 
This configuration is a strong evidence in favour of the trig- 



gering star formation process descripted by the "collect and 
collapse" model. In all of the preceding examples, we have 
an Hn region generated by an O-type star or a cluster of 
massive stars. 

The region around Shell B presents a similar interstel- 
lar scenario as the preceding examples, e.i an optical shell 
bounded by a strong arc-like feature of PAH emission which 
are surrounding by molecular material (Clouds A, E, and 
F). In such scenario, star formation can be induced by the 
action of the massive stars through the ionizing flux and stel- 
lar winds. Indeed, the presence of YSOs towards the molec- 
ular clouds linked to Sh2-132 indicates that star formation 
is very active in this region of the Galaxy (Vasquez et al. in 
preparation) . 
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